Finite driving rate and anisotropy effects in landslide modeling 
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In order to characterize landslide frequency-size distributions and individuate hazard scenarios 
and their possible precursors, we investigate a cellular automaton where the effects of a finite driving 
rate and the anisotropy are taken into account. The model is able to reproduce observed features of 
landslide events, such as power-law distributions, as experimentally reported. We analyze the key 
role of the driving rate and show that, as it is increased, a crossover from power-law to non power-law 
behaviors occurs. Finally, a systematic investigation of the model on varying its anisotropy factors 
is performed and the full diagram of its dynamical behaviors is presented. 
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I. INTRODUCTION 

Over the last two decades, the evidence of power laws 
in frequency-size distributions of several natural hazards 
such as earthquakes , volcanic eruptions , forest fires 
0- LH HI and landslides @, has suggested a relationship 
between these complex phenomena and self-organized 
criticality (SOC) 0. The idea of SOC 0, applied to 
many media exhibiting avalanche dynamics |3, ll(J | , refers 
to the tendency of natural systems to self-organize into a 
critical state where the distribution of event sizes is rep- 
resented by a power law with an exponent a, which is 
universal in the sense that it is robust to minor changes 
in the system. Generally, the nature of a critical state is 
evidenced by the fact that the size of a disturbance to the 
system is a poor predictor of the system response. Let 
us consider storms as perturbations for natural slopes. 
Large storms can produce large avalanches, but also small 
storms sometimes can do it. On the other hand, small 
storms usually do not produce any avalanche, but also 
large storms may not cause any avalanching phenomena. 
Moreover, avalanches triggered by small storms can be 
larger than those triggered by large storms. The un- 
predictability of the sizes of such system responses to 
incremental perturbations and the observed power-law 
statistics could be the exhibition of self-organized criti- 
cal behavior in most natural avalanches. However, the 
idea of understanding power-law distributions within the 
framework of SOC is not the only one. Recently, in or- 
der to reproduce the range of the power-law exponents 
observed for landslides, some authors have introduced a 
two-threshold cellular automaton, which relates landslide 
dynamics to the vici nity of a breakdown point rather than 
to self-organization pd| . 

In this paper, we report an accurate investigation of a 
cellular automaton model, which we have recently pro- 
posed to describe landslide events and, specifically, their 
frequency-size distributions |l2j |. In particular, we dis- 
cuss the role of a finite driving rate and the anisotropy ef- 
fects in our non-conservative system. It has been pointed 
out by several authors that the driving rate is a parame- 
ter that has to be fine tuned to zero in order to observe 
criticality O 13, E^- We notice that the limit of zero 



driving rate is only attainable in an ideal chain reaction, 
therefore finite rates of external drives are essential ingre- 
dients in the analysis of the dynamics of real avalanche 
processes. We show that increasing the driving rate 
the frequency-size distribution of landslide events evolves 
continuously from a power-law (small driving rates) to 
an exponential (Gaussian) function (large driving rates). 
Interestingly, a crossover regime characterized by a max- 
imum of the distribution at small sizes and a power-law 
decay at medium and large sizes is found in the interme- 
diate range of values of the driving rate for a wide range 
of level of conservation. Power-law behaviors are robust 
even though their exponents depend on the system pa- 
rameters (e.g., driving rate and level of conservation, see 
below) . 

Although the critical nature of landslides is not fully 
assessed and many authors believe that deviations from 
power la w ap pear to be systematic for small landslides 
data 0, Il7| . results from several regional landslide in- 
ventories show robust power-law distributions of medium 
and large events with a critical exponent a ~ 2.5 ± 0.5 
. The variation in the exponents of landslide size dis- 
tributions is larger than in the other natural hazards that 
exhibit scale-invariant size statistics [lg. Whether this 
variation of a is caused by scatter in the data or because 
different exponents are associated with different geology, 
is an important open question, which we may contribute 
to address. 

The model we analyze describes the evolution of a 
space and time dependent factor of safety field. The 
factor of safety (FS) is defined as the ratio between re- 
sisting forces and driving forces. It is a complicate func- 
tion of many dynamical variables (pore water pressure, 
lithostatic stress, cohesion coefficients, etc.) whose rate 
of change is crucial in the characterization of landslide 
events. A landslide event may include a single landslide 
or many thousands. We investigate frequency-size dis- 
tributions of landslide events by varying the driving rate 
of the factor of safety. Although our probability den- 
sity distributions are lacking of a direct comparison with 
frequency-size distributions of real landslides they repro- 
duce power-law scaling with an exponent very close to 
the observed values. Moreover, they allow us to get in- 
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sight into the difficult problem of the determination of 
possible precursors of future events. 

The paper is organized as follows. In the next Section, 
we present the model and briefly discuss the differences 
between our approach and previous cellular automata 
models that have been recently introduced to character- 
ize landslide frequency-size distributions. In Section III, 
we report numerical results obtained by a systematic in- 
vestigation of the effects of a finite driving rate on the 
frequency-size distribution. The values of the exponent 
of the power-law decay are given as a function of the 
driving rate and the level of conservation. An accurate 
analysis of the spatial distribution of the values of the fac- 
tor of safety by varying the driving rate provides useful 
information for quantifying hazard scenarios of possible 
avalanche events. In Section IV, we analyze the role of 
anisotropic transfer coefficients, which control the prop- 
agation of the instability. We summarize our results in 
a phase diagram that shows the location of power-law 
and non power-law scaling regions in the anisotropy pa- 
rameter space. Conclusions are summarized in Section 
V. 



II. THE MODEL 

The instability in clays often starts from a small re gion , 
destabilizes the neighborhood and then propagates [19j. 
Such a progressive slope failure recalls the spreading of 
avalanches in the fundamental models of SOC. The term 
self-organized criticality (SOC) was coined by Bak, Tang 
and Wicscnfcld to describe the phenomenon observed in a 
particular cellular automaton model, nowadays known as 
the sandpile model [8j . In the original sandpile model, the 
system is perturbed externally by a random addition of 
sand grains. Once the slope between two contiguous cells 
has reached a threshold value, a fixed amount of sand is 
transferred to its neighbors generating a chain reaction 
or avalanche. The non-cumulative number of avalanches 
Na with area A satisfies a power-law distribution with 
a critical exponent a = 1 20], which is much smaller 
than the values of the power-law exponents observed for 
landslides [UEl- Few years later the paper of Bak et al. 

Olami, Feder and Christensen (OFC) recognized the 
dynamics of earthquakes as a physical realization of self- 
organized criticality and introduced a cellular automaton 
that gives a good prediction of the Gutenberg- Richter law 
. Such a model, whose physical background belongs to 
the Burridge-Knopoff spring-block model |2l|, is based 
on a continuous dynamical variable which increases uni- 
formly through time till reaches a given threshold and re- 
laxes. This means that the dynamical variable decreases, 
while a part of the loss is transferred to the nearest neigh- 
bors. If this transfer causes one of the neighbors to reach 
the threshold value, it relaxes too, resulting in a chain 
reaction. OFC recognized that the model still exhibits 
power-law scaling in the non-conservative regime, even if 
the power-law exponent strongly depends on the level of 



conservation. 

In this paper, we investigate the role of a finite driv- 
ing rate and of anisotropy in a non-conservative cellular 
automaton modeling landslides In such a model, we 
sketch a natural slope by using a square grid where each 
site i is characterized by a local value of the safety factor 
FSi. In slope stability analysis, the factor of safety, FS, 
against slip is defined in terms of the ratio of the maxi- 
mum shear strength r max to the disturbing shear stress 

T 

FS =T^. (l) 

The limited amount of stress that a site can support is 
given by the empirical Mohr-Coulomb failure criterion: 
r max = c + (er — u) tan <f>, where a is the total normal 
stress, u is the pore-fluid pressure, cf> is the angle of in- 
ternal friction of the soil and c is the cohesional (non- 
frictional) component of the soil strength • If FS > 1 , 
resisting forces exceed driving forces and the slope re- 
mains stable. Slope failure starts when FS = 1. Since 
a natural slope is a complex non-homogeneous system 
characterized by the presence of composite diffusion, dis- 
sipative and driving mechanisms acting in the soil (such 
as those on the water content), we consider time and site 
dependent safety factor FSi and treat the local inverse 
factor of safety = 1 /FSi as the non-conserved dynam- 
ical variable of our cellular automata model 0] . 

The long-term driving of the OFC model is, here, re- 
placed by a dynamical rule which causes the increases 
of ei through the time with a finite driving rate v: 
ei(t + At) = ei(t) + i/At. Such a rule allows us to simulate 
the effect on the factor of safety of different complex pro- 
cesses which can change the state of stress of a cell. The 
model is driven as long as < 1 on all sites i. Then, 
when a site, say i, becomes unstable (i.e., exceeds the 
threshold, e t h = 1) it relaxes with its neighbors accord- 
ing to the rule: 

&i * 0, 6™ ^ Grm ~t~ fnn^ii (2) 

where nn denotes the nearest neighbors of site i and /„„ 
is the fraction of ei toppling on nn. This relaxation rule 
is considered to be instantaneous compared to the time 
scale of the overall drive and lasts until all sites remain 
below the threshold. When reaches the threshold value 
1 and relaxes, the fraction of moving from the site i 
to its "downward" (resp. "upward") neighbor on the 
square grid is fd (resp. f u ), as /; = /,■ is the fraction 
to each of its "left" and "right" neighbors. The trans- 
fer parameters /„„ are chosen in order to individuate a 
privileged transfer direction: we assume /„ < fd and 
fl = fr < fd- We notice that the model reproduces fea- 
tures of the OFC model for earthquakes in the limit case 
v = and /„„ = / < 0.25. A detailed analysis of the 
model on varying the transfer coefficients /„„ is reported 
in Sec. IV. 

Since many complex dissipative phenomena (such 
as evaporation mechanism, volume contractions, etc. 
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(23|) contribute to a dissipative stress transfer in 
gravity-driven failures, we study the model in the non- 
conservative case C — J2nn fnn < lj which makes our 
approach different from previous ones within the frame- 
work of SOC [lg. The conservation level, C, and the 
anisotropy factors, which we consider here to be uni- 
form, are actually related to local soil properties (e.g., 
lithostatic, frictional and cohesional properties), as well 
as to the local geometry of the slope (e.g., its morphol- 
ogy) . The rate of change of the inverse factor of safety, 
v, induced by the external drive (e.g., rainfall), in turn 
related to soil and slope properties, quantifies how the 
triggering mechanisms affect the time derivative of the 
FS field. 

Recently, in order to reproduce the range of the power- 
law exponents observed for landslides, several authors 
have used two-threshold cellular automata, which re- 
late landslide dynamics to self-organization [2J| or to the 
vicinity of a breakdown point |ll| . In the first approach 
|24) . a time-dependent criterion for stability, with a not 
easy interpretation in terms of governing physics, pro- 
vides a power-law exponent close to 2 without any tuning. 
Therefore, this approach does not explain the observed 
variability of a. In Ref. the range of a is found 

by tuning the maximum value of the ratio between the 
thresholds of two failure modes, the shear failure and 
the slab failure. However, the frequency-size distribu- 
tion of avalanches is obtained by counting only clusters 
where shear failures have occurred, considering conserva- 
tive transfer processes between adjacent cells with a dif- 
ferent number of nearest neighbors. In this paper, the in- 
vestigation of our non-conservative cellular automaton is 
mainly devoted to the characterization of landslide event 
dynamics on varying the driving rate in order to analyze 
different hazard scenarios. 



III. THE EFFECT OF A FINITE DRIVE ON 
FREQUENCY-SIZE DISTRIBUTIONS 
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FIG. 1: Noncumulative frequency-size distributions 011 a 
64 x 64 grid corresponding to four values of the driving rate. 
The logarithm of the normalized number of model events, 
logia[P(s)}, in which a specified number of different cells, s, 
become unstable, is given as a function of logio(s). We show 
the case C = 0.4, f u /f d = 2/3 and f t /f d = 5/6. 
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Frequency-size distributions give the number of land- 
slides (events) as a function of their size. In FigQJwe show 
the non-cumulative frequency-size distributions obtained 
for different values of the driving rate in the anisotropic 
non conservative case C = 0.4, with f u /fd = 2/3 and 
fl/fd = 5/6. The curves are obtained for a square lat- 
tice of size 64 x 64. We considered both cylindrical (open 
along the vertical axis and periodical along the horizontal 
axis) and open boundary conditions, which we checked 
differ in the slopes of the distribution curves for less than 
1%. 

In the limit of vanishing driving rate, the distribution 
of events, P(s), is similar to that of the two-dimensional 
isotropic OFC model for a fixed value of the level of con- 
servation: a power law characterized by a critical expo- 
nent a, P{s) ~ s~ Q , followed by a system finite-size de- 
pendent exponential cutoff Q . As discussed in Ref. 0] , 
by increasing the driving rate u, the probability distri- 



FIG. 2: Noncumulative frequency-size distributions and cor- 
responding Gaussian distributions on a 64 x 64 grid for two 
values of the driving rate. Squares and triangles show the nor- 
malized number of model events, P(s), in which a specified 
number of different cells, s, become unstable as a function of 
s. The dashed and the solid lines are Gaussian distributions 
obtained for the same mean value and the standard devia- 
tion of the frequency-size distributions. We consider the case 
C = 0.4, f u /f d = 2/3 and fi/f d = 5/6. 



bution develops a maximum, which shifts towards larger 
events with v. On the left side of the maximum, a power- 
law decay, with exponent j3 (see the Fig. QJ, seems to ap- 
pear for small landslide sizes. However, the few available 
data do not allow to distinguish log-log linear shape and 
an exponential one 17]. On the right side of the maxi- 
mum of the distribution, the power-law regime remains 
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FIG. 3: Top: Positive power-law exponent j3 as a function of 
the driving rate v. Bottom: Negative power-law exponent a 
as a function of v. The values of the exponents are obtained 
for C = 0.4, U/U = 2/3 and fi/f d = 5/6. 
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FIG. 4: Top: Positive power-law exponent j3 as a function 
of the level of conservation C. Bottom: Negative power-law 
exponent a as a function of C. The values of the exponents 
are obtained for v = 0.003, fu/fd = 2/3 and fi/fd = 5/6. 



until, by increasing v 1 the distribution continuously mod- 
ifies in a bell-shaped curve. Fig|2 shows the crossover of 
the probability distribution from power-law to Gaussian 
on increasing the driving rate. 

The behavior of the power-law exponents on varying v 
is shown in FigEJ Interestingly, the values of the power- 
law exponent a are very close to those experimentally 
reported As one can see, a increases until the value 
of the driving rate sensitively modifies the shape of the 
distribution with the appearance of a maximum. We 
find that the regime with a non-monotonic frequency-size 
distribution is robust to changes in system parameters. 
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FIG. 5: Size, s, of equiprobable events (i.e., corresponding 
to the same probability of occurrence) as a function of the 
driving rate. Top: Sizes of model events with P(s) = 10~ 2 . 
Bottom: Sizes of model events with P(s) = 10 -6 . In the 
top panel, only the largest sizes of equiprobable events are 
plotted. 



In particular, it can be found for v E [10~ 4 , 1CP 2 ] in the 
whole range C € [0.4, 0.8]. 

The values of the exponents a and f3 sensitively change 
by varying another important parameter of the model 
that is the level of conservation C, which represents the 
non-conservative redistribution of the load of failing cells. 
The effect of C at finite driving rate is comparable to that 
obtained at v = |lj, as shown in FigQ] (see also Ref. 

G2)- 

Let us come back to Fig^ in order to highlight the 
role of the driving rate in the characterization of land- 
slide events. From the distributions of Fig 2] we see that 
large events can have comparatively high probabilities in 
the power-law regime (i.e., at small driving rates v) with 
respect to the Gaussian regime (i.e., larger v). Such a fea- 
ture could allow to reproduce the observation that some- 
times mild rainfalls produce landslides as large as those 
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FIG. 6: (top row) Snapshots of a landslide event of size s = 230 on a 64 x 64 grid, for four values of the driving rate (from left 
to right): v = 10~ 5 , v — 1CP 3 , v = 5 ■ 10~ 3 and v = 5 • 1CF 2 . The 230 black cells are those that have reached the instability 
threshold. The simulations are done in the case C = 0.4, fu/fd = 2/3 and fi/fd = 5/6. (bottom row) Snapshots of the factor 
of safety corresponding to the stable configurations reached after the avalanches shown in top row. The values of the factor of 
safety have been associated to ten levels of a gray scale from white to black, in order to measure the distance of a cell from its 
instability condition: the darker the color, the farther is the cell from the instability threshold. 



triggered by intense rainfalls. To underline the nonlinear 
behavior of the system with v, in Fig|3]we plot the sizes 
of avalanches with the same probability of occurrence for 
different values of the driving rate. Such sizes are the in- 
tersections of horizontal straight lines with the distribu- 
tion curves of Fig^ When we consider small events (i.e., 
P(s) as large as 10 -2 ), we find that the size of equiprob- 
able events increases with the driving rate. Thus, for 
equiprobable events with high probability of occurrence 
the system response is essentially linear with the driv- 
ing rate. Instead, for large events (i.e., P(s) as small as 
10~ 6 ), the size of equiprobable events has a maximum as 
a function of v. Thus, it appears that, for a given range 
of v values, the size of events caused by a slow rate of 
changes of the factor of safety can be larger than the size 
of avalanches triggered by a faster rate. Moreover, an 
evidence is found for the existence of the most dangerous 
value of v, for which the size of the system response has 
a maximum. 



A. Hazards scenarios 

The detection of possible precursors of a landslide 
event is a crucial step to achieve hazard reduction. In 
order to get insight into this difficult problem, we visual- 
ize the structure of a typical landslide event for different 
values of the driving rate v. In the top row of Fig0 
we report on the 64 x 64 grid a typical avalanche of size 
s = 230 on increasing the value of v from left to right. 
The 230 black cells are those that have reached the in- 
stability threshold. As we can see, compact landslides 
are the characteristic response of a system governed by 
power-law statistics as it happens at small v (also when 
a maximum in the frequency-size distribution develops). 
Such a response is typical of systems with SOC behav- 
ior [2(|. By increasing the driving rate, compact clusters 
survive until power-law regime disappears. As the sys- 
tem enters the non power-law regime, the relevance of 
domino effects drastically drops and landslide events are 
characterized by many tiny independent clusters. 

In the bottom row of FiglHI we show the distribution of 
the factor of safety FSi = 1/e,; for the cases correspond- 
ing to the upper panels |2q . The distribution on the 
grid of the values shows how the spatial correlations in 
the system crucially affect landslide structures. In order 
to measure the distance of a cell from its instability con- 
dition and to visualize the correlated areas (regions with 
similar values of FSi), the values of FSi have been asso- 
ciated to ten levels of a gray scale from white to black: 



the darker the color, the farther is the cell from the insta- 
bility threshold. In the snapshots of FigEJ it is possible 
to recognize as dark areas the avalanches shown in the 
corresponding upper grids. In particular, the dark areas 
typically are related to previous landslide events, whereas 
the lighter areas indicate regions of future events. We 
notice that in the power law regime (i.e. small v) even 
a very small perturbation (say, a drop of water) at one 
single point can trigger huge system responses. Instead, 
in the non power-law regime (i.e. large v) large-scale 
correlations are absent; here large events trivially occur 
just because the strong external driving rate makes likely 
that many cells simultaneously approach the instability 
threshold. Thus, the detection of patterns of correlated 
domains in investigated areas results to be a crucial tool 
to identify the response of the system to perturbations, 
i.e., to hazard assessment. 

It is worth noticing that the average value of the fac- 
tor of safety on the grid cells < FS > and its fluctuation 
< AFS 2 > are very similar in the four cases showed in 
FigEl encompassing a broad spectrum of v values. Inter- 
estingly, the probability distribution of FS on the grid 
sites, independently of the driving rate, is well approxi- 
mated by a Gaussian distribution. This suggests that a 
measure of just an average safety factor on the investi- 
gated area could provide only a very partial information 
about the statistics governing the considered landslide 
events. 



IV. ANISOTROPY EFFECTS ON 
FREQUENCY-SIZE DISTRIBUTIONS 

In the previous sections, we have investigated the prop- 
erties of the model on varying the driving rate v at fixed 
values of the anisotropic ratios f u /fd and fi/fd = fr/ fd- 
Such parameters control how the instability propagates 
downward and, therefore, they are complicate functions 
of the topography and geology of a specific area. In this 
Section, we are interested in the analysis of the model 
when such transfer coefficients vary. 

It is well-known that the one-dimensional version of 
the sandpile and OFC models are characterized by non 
power-law scaling 0. Thus, we expect that, for small val- 
ues of the ratio fi/ fd, the frequency-size distribution docs 
not show power-law behavior. Viceversa, for f u / fd ~ 1 
and fi/fd ~ 1) we expect to have power-law distribu- 
tions, as in the OFC model, which corresponds in our 
model to the limit / = f u = fi = f r = fd and v = 0. 

The diagram of Fig [7] summarizes the different regimes 
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FIG. 7: Phase diagram in the f u /fd vs. fi/fd plane. The 
lines divide the parameter space in three regions: a power- 
law scaling region (PL) , a non power- law scaling region (NPL) 
and a finite-size oscillation region (FSO) where peaks of the 
distribution commensurate with the size of the grid appear. 
The results are obtained for v = 10~ 4 and C = 0.4. 
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FIG. 8: Noncumulative frequency-size distributions on a 64 x 
64 grid for f u /fd = and fi/fd = 0.3 (triangles) and fi/f d = 
0.55 (squares). The inset shows the distribution curves in a 
log-linear scale for the range of values where commensurate 
peaks are observed. The vertical lines mark multiples of the 
system size. The results are obtained for v = 10~ 4 and C = 
0.4. 



found in our simulations at fixed values of the driving 
rate v — 10~ 4 and the level of conservation C = 0.4. We 
find that, on varying the anisotropic ratios, the parame- 
ter space is divided in three regions: i) a power-law region 
(PL) characterized by power-law frequency-size distribu- 
tions for large values of the anisotropic ratios, it) a non 
power-law region (NPL) for the whole range of values of 
the anisotropic ratio f u / fd (which controls the redistri- 
bution of load in the vertical direction) and small values 
°f fi/fd, Hi) a finite-size oscillation region (FSO) where 
the frequency-size distribution is characterized by peri- 
odic peaks, which appear for integer multiples of the grid 
size L = 64. 

As expected, we find power-law and non power-law 
behaviors for large and small values of the anisotropic 
ratios, respectively. In particular, we find that the value 
of the critical exponent a of the power-law distribution 
slightly increases with fu/fd, for a fixed value of fi/fd and 
decreases with fi/fd for a fixed value of f u /fd- However, 
the changes in a are negligible. 

It is worth noticing that, even if the ratio fij fd is quite 
large, we find that the frequency-size distribution does 
not show a power law decay when the value of f u / fd is 
small, (see Fig[7|). Indeed, the probability distribution 
develops a finite number of peaks which are commensu- 
rate with the size of the grid. FiglHl shows the frequency- 
size distribution for two different values of the anisotropic 
ratios in the FSO region of the phase diagram. By in- 
creasing fi/ fd, the peaks of the distributions turn down 
as long as they disappear. We m that several peaks in the 
distribution of avalanche sizes are obtained in Ref. [2^| . 
The authors vary the convexity of the driving 7 and the 
level of conservation e in the isotropic case: C = I when 
e = 0.25. They find an intermediate region of the phase- 
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FIG. 9: Noncumulative frequency-size distributions on a 64 x 
64 grid for fu/fd = 1 and fi/fd = 0.3 (open dots) and fi/f d = 
0.7 (squares). The results are obtained for v = 10 -4 and 
C = 0.4. 

diagram e vs. 7 where the frequency-size distribution is 
characterized by peaks that scale with different powers of 
the system size L. Fixing the convexity of the driving 7 
(7 = for a uniform driving) , for large values of the level 
of conservation, C > 0.6, the peaks disappear and they 
get a power-law decay with an exponential cutoff. Con- 
versely, we find commensurate peaks in the FSO region 
in the whole range C 6 [0.4,0.8]. 

We notice that an analysis of the anisotropic case of 
the OFC model is made in Ref.[27j where the authors 
introduce only two transfer coefficients ot\ = fi = f r and 
a 2 = fu = fd and control the degree of anisotropy by 
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changing the ratio a,\/a.2, while keeping the level of con- 
servation constant. They find that the anisotropy has 
almost no effect on the power-law exponent while the 
scaling exponent, expressing how the finite-size cutoff 
scales with the system size, changes continuously from 
a two-dimensional to a one dimensional scaling of the 
avalanches |27|. Varying the anisotropic ratio oti/a.2 in 
the range [0, 1] is equivalent to consider the straight line 
fu/fd = 1 in the phase diagram of Fig0 As in Ref. |27| . 
we find that on moving along the line f u /fd — 1, the 
changes in the power-law exponent are negligible. How- 
ever, differently Ref. [27]], we find a crossover in the 
frequency-size distribution behavior from power-law to 
non power-law (see Fig. EJ. We attribute such a differ- 
ent result to the finite driving rate. 

In conclusion, our analysis shows that only a finite 
range of values of the anisotropic transfer coefficients 
can supply power-law distributions. This characteriza- 
tion provides insight into the difficult determination of 
the complex and non-linear transfer processes that occur 
in a landslide event. 



tomaton based on a dissipative factor of safety field. We 
have found that the value of the driving rate, which de- 
scribes the variation rate of the factor of safety due to 
external perturbations, has a crucial role in determining 
landslide statistics. In particular, we have shown that, as 
the driving rate increases, the frequency-size distribution 
continuously changes from power-law to gaussian shapes, 
offering the possibility to explain the observed rollover of 
the data for small landslides. The values of the calcu- 
lated power-law exponents are in good agreement with 
the observed values. Moreover, the analysis of the model 
on varying the driving rate suggests the determination 
of correlated spatial domains of the factor of safety as 
a useful tool to quantify the severity of future landslide 
events. 

As concerns the effects of anisotropic transfer coef- 
ficients, which control the non-conservative redistribu- 
tion of the load of failing cells, we have found that the 
power-law behavior of the frequency-size distribution is 
a feature of the model only in a limited region of the 
anisotropy parameter space. 



V. CONCLUDING REMARKS 

Explanation of the power-law statistics for landslides is 
a major challenge, both from a theoretical point of view 
as well as for hazard assessment. In order to characterize 
frequency-size distributions of landslide events, we have 
investigated a continuously driven anisotropic cellular au- 



Acknowledgments 

E. Piegari wishes to thank A. Avella for stimulating 
discussions and a very friendly collaboration. This work 
was supported by MIUR-PRIN 2002/FIRB 2002, SAM, 
CRdC-AMRA, INFM-PCI, EU MRTN-CT-2003-504712. 



[1] Z. Olami, H. J. S. Feder, K. Christensen, Phys. Rev. Lett. 

68, 1244 (1992). 
[2] T. Simkin, Annu. Rev. Earth Planet. Sci. 21, 427 (1993). 
[3] P. Bak, K. Chen, C. Tang, Phys. Lett. A. 147, 297 

(1992). 

[4] R. Pastor-Satorras, and A. Vespignani, Phys. Rev. E 61, 
4854 (2000). 

[5] D. L. Turcotte, B. D. Malamud, F. Guzzetti, P. Reichen- 
bach, Proc. Natl. Acad. Sci. U.S.A. 99, 2530 (2002). 

[6] C. Dussauge, JR. Grasso, and A Helmstetter, J. Geo- 
phys. Res. 108, 2286 (2003). 

[7] P. Bak, How Nature Works - The Science of Self- 
Organized Criticality, (Copernicus, Springer- Verlag, New 
York, 1996). 

[8] P. Bak, C. Tang and K. Wiesenfeld, Phys. Rev. Lett. 59, 

381 (1987); Phys. Rev. A 38, 364 (1988). 
[9] H. J. Jensen, "Self-Organized Criticality: emergent com- 
plex behavior in physical and biological systems" (Cam- 
bridge University Press, Cambridge, 1998). 

[10] D. L. Turcotte, Rep. Prog. Phys. 62, 1377 (1999). 

[11] J. Failletaz, F. Louchet, JR. Grasso, Phys. Rev. Lett. 
93, 208001 (2004). 

[12] E. Piegari, V. Cataudella, R. Di Maio, L. Milano, M. 
Nicodemi, Geophys. Res. Lett, in press. 

[13] D. Hamon, M. Nicodemi and H.J. Jensen, Astron- 
omy^ Astrophysics 387, 326 (2002). 

[14] R. Dickman, A. Vespignani, S. Zapperi, Phys. Rev. E 57, 



5095 (1998). 

[15] D. Sornette, Critical Phenomena in Natural Sciences, 
Chaos, Fractals, Self- organization and Disorder: Con- 
cepts and Tools, (Springer Series in Synergetics, Heidel- 
berg, 2004). 

[16] F. Brardinoni, M. Church, Earth Surf. Process. Land- 
forms 29, 115 (2004). 

[17] B. D. Malamud, D. L. Turcotte, F. Guzzetti, P. Reichen- 
bach, Earth Surf. Process. Landforms 29, 687 (2004). 

[18] S. Hergarten, Natural Hazards and Earth System Sci- 
ences 3, 505 (2003). 

[19] L. Bjerrum, J. Soil Mech. Fdns. Div. Am. Soc. Civ. En- 
gnrs. 93, 3 1967. 

[20] L. P. Kadanoff, S. R. Nagel, L. Wu, S. M. Zhou, Phys. 
Rev. A 39, 6524 (1989). 

[21] R. Burridge and L. Knopoff, Bull. Seismol. Soc. Am. 57, 
341 (1967). 

[22] K. Terzaghi, Geothecnique 12, 251 (1962). 

[23] D. G Fredlund, H. Rahardjo, "Soil Mechanics for Unsa- 

tured Soils" (Wiley-Interscience, New York, 1993). 
[24] S. Hertgarten and H. J. Neugebauer, Phys. Rev. E 61, 

2382 (2000). 

[25] In order to visualize the FSi distribution, we introduced 
a lower threshold e m i„ = 10 -3 and checked that the re- 
sults do not change. 

[26] L. Pietronero, W. R. Schneider, Phys. Rev. Lett. 66, 2336 
(1991). 



8 

[27] K. Christensen and Z. Olami, Phys. Rev. A 46, 1829 Rev. Lett. 74, 118 (1995) 

(1992) 

[28] A. Corral, C. J. Perez, A. Diaz-Guilera, A. Arenas, Phys. 



This figure "fig6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0603100vl 



